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Abstract 

The three- leg ladder has one odd-parity and two even-parity channels. At 
low doping these behave quite differently. Numerical calculations for a t — J 
model show that the initial phase upon hole doping has two components - 
a conducting Luttinger liquid in the odd-parity channel, coexisting with an 
insulating (i.e. undoped) spin liquid phase in the even-parity channels. This 
phase has a partially truncated Fermi surface and violates the Luttinger the- 
orem. This coexistence of conducting fermionic and insulating paired bosonic 
degrees of freedom is similar to the recent proposal of Geshkenbein, Ioffe, and 
Larkin (Phys. Rev. B55, 3173 (1997)) for the underdoped spin-gap normal 
phase of the cuprates. A mean field approximation is derived which has many 
similarities to the numerical results. One difference however is an induced 
hole pairing in the odd-parity channel at arbitrary small dopings, similar to 
that proposed by Geshkenbein, Ioffe, and Larkin for the two-dimensional case. 
At higher dopings, we propose that a quantum phase transition will occur as 
holes enter the even-parity channels, resulting in a Luther-Emery liquid with 
hole pairing with essentially d-wave character. In the mean field approxima- 
tion a crossover occurs which we interpret as a reflection of this quantum 
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phase transition deduced from the numerical results. 
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I. INTRODUCTION 



The properties of electrons confined to ladders with various numbers of legs have been 
investigated by many groups in the past few years |1J . The reason for this interest lies both in 
the unusual properties of the ladder systems and the possibility to realize ladder structures 
in the cuprates, but also in the insight these systems give to the full two-dimensional problem 
of a square lattice. Two different approaches have been taken. One is based on the weak- 
coupling limit and uses renormalization group methods to analyze the multi-leg ladders 
. A very complete analysis of this type for N-leg ladders in the Hubbard model has 
recently been given by Lin, Balents and Fisher |7j]. 

A second approach is more numerical and examines the strong-coupling limit described 
mostly by the t-J model. Recent progress on the loop algorithm for Monte Carlo calculations 
has allowed large undoped systems described by a Heisenberg model to be investigated down 
to very low temperatures. However, when doped holes are introduced, the fermion sign 
problem prevents one from using this method, and other methods must be employed, e.g. 
using Lanczos techniques to diagonalize relatively small systems |||| or the new density 
matrix renormalization group method (DMRG) to obtain the groundstate of large systems 

In this paper we examine the case of the lightly doped t — J three-leg ladder. This case 
is specially interesting because in a certain sense it combines the contrasting properties of 
a single chain [|ll|!0] and a two-leg ladder |§. These behave quite differently, both when 
undoped or when lightly doped, so it is of great interest to follow the evolution of the three-leg 
ladder in this regime. In particular, the evolution of the Fermi surface as the Mott insulating 
state is approached is very different in the different transverse channels. As we shall discuss 
further below, this leads to a region where the Fermi surface is truncated in two channels, 
but remains in one channel - a behavior which clearly violates the Luttinger theorem. We 
shall argue that this presence of transverse channels which behave quite differently helps 
us to make inferences for the limit of the two-dimensional plane which of course can be 
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represented as the limit of very many channels or patches on the Fermi surface. The lightly 
doped limit of the three-leg ladder can serve as a simpler analog for the underdoped spin 
gap region of the cuprates, which has attracted so much attention, if one assumes that the 
patches or channels near (±7T, 0) and (0, ±7r) are truncated and show a spin gap, while those 
near (±7r/2, ±7r/2) are gapless ]T5| . 

The t-J three- leg ladder Hamiltonian is given by 

H= -tJ2j2P(c}^c J+1 ^ + H.c.)P 
-t'Y.jZP{c)^c hU+ ^ + H.c.)P 

j,cr u=l 

3 1 

j v=\ 



j u=l 4 



(1) 



where j runs over L rungs, cr(=t, |), and v are spin and leg indices. The t-J three-leg ladder 
is sketched in Fig. 1. The first two terms are the kinetic energy (P is a projection operator 
which prohibits double occupancy), and the last two exchange couplings J (J') act along 
the legs (rungs). Periodic or antiperiodic boundary conditions (PBC, APBC) are used along 
the legs. 
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FIG. 1. The t-J ladder with three legs and L rungs. The couplings along the legs are t and J, 
and those along the rungs t' and J' . 



The paper is organized as follows. In the next section we recapitulate briefly the known 
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results in the undoped limit described by a Heisenberg Hamiltonian. Then in section 3 we 
discuss the case of a single doped hole, using the results of a Lanczos diagonalization and 
also the earlier DMRG results by White and Scalapino ||10|| . The low-energy properties are 
described by a single Luttinger liquid channel in contact with an insulating spin liquid (ISL). 
Next in section 4 we consider states with two and more holes. In the presence of a finite hole 
density, there are two possibilities - either all the holes repel each other and enter the single 
Luttinger liquid channel, or at some density, the other channels are also populated with 
holes. Since doping a resonating valence bond (RVB) spin liquid leads to a Luther-Emery 
liquid ||, this will cause a qualitative change in the physical properties. There will be a 
critical hole density, 5 C , which controls the transition between the low-density phase with 
only a Luttinger liquid, in contact with a spin liquid, to the case with both Luttinger and 
Luther-Emery liquids. In section 5 we develop a mean field approach to the t-J three-leg 
ladder. In this mean field approximation, as we shall see, the critical density, 8 C , is not finite 
but arbitrarily small, whereas the numerical results give a finite 5 C for values of J/t ~ 0.5. 
In the final section we draw some conclusions and discuss the relationship to the planar 
two-dimensional case. 



II. UNDOPED CASE: THE HEISENBERG THREE-LEG LADDER 

In this section we recapitulate briefly the known results for three-leg ladders in the 
Heisenberg model. In this case the fermion sign problem does not occur, and very accurate 



quantum Monte Carlo calculations have been carried out by Frischmuth et al. JL4], and by 



Greven et al. [|T^]. The low-energy properties can be mapped onto an effective single S=l/2 



antiferromagnetic (AF) Heisenberg chain model, although the starting model has three spins 
per set of ladder rungs. This arises because the additional spin degrees of freedom have a 
spin gap, and enter with an excitation energy ~ J'. The spinon velocity of the effective 
model in fact is hardly changed from the value of the nearest neighbor Heisenberg chain 
with exchange constant J. However the energy scale parameter which controls the spinon- 
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spinon interactions and the size of logarithmic corrections to the uniform spin susceptibility 



at finite temperatures is greatly changed. Frischmuth et al. [OJ interpret this result in terms 
of a Heisenberg model for the effective chain with longer range and unfrustrated effective 
interactions which act to enhance the AF correlations. 

The key point however is the fact that the low-energy Hilbert subspace is greatly reduced 
so that not three chains but only a single chain model is active at low energies, or in other 
words only a single S=l/2 degree of freedom per rung of three spins remains in the low- 
energy region. The remaining spin degrees of freedom are gapped in the same way as the 
spin liquid of a two-leg ladder. 
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FIG. 2. (a) and (b) Spinon excitation spectra for the isotropic (3 x 7)- and (3 x 8)-site Heisenberg 
ladders. The circles correspond to the even parity and the squares to the odd parity channels 
respectively, (a) Single spinon excitation spectrum (S = 1/2) in a (3 x 7)-cluster. (b) Two-spinon 
excitation spectrum (S = 1) in a (3 x 8)-cluster. (c) Spinon-velocity, v s , in (3xL)-Heisenberg 
ladders. v s extrapolates with high accuracy to the value 1.81J (diamond), obtained by QMC 
calculations [13]. 



This behavior is clearly observed in the numerical Lanczos diagonalization of (3 x 7)- and 
(3 x 8)-clusters, shown in Fig. 2. As Bares et al. [12| emphasized in their study of the exactly 



solvable supersymmetric t-J model in one dimension, the spinon and holon dispersions can 
be directly obtained in the limit of vanishing hole density by examining the groundstate 
of an appropriately chosen finite chain. Thus if one takes the case of one undoped chain 
with an odd number of sites and PBC, then the groundstate for each total wavevector 
along the chain gives the dispersion of a single spinon. This groundstate manifold has a 
total spin quantum number S=l/2. For the present case of a three- leg ladder we show the 
dispersion for a (3 x 7)-sample with PBC in Fig. 2(a). The single spinon with S=l/2 has 
odd parity with respect to reflection about the center leg and disperses with a bandwidth 
~ | J and minima at ±|. The corresponding one-spinon excitation spectrum of a 7-site 
single Heisenberg chain is shown in Fig. 3(a). To estimate finite size effects, the calculation 
was repeated for a 19-site Heisenberg chain (Fig. 3(c)). We observe that the qualitative 
features, i.e. the positions of extrema and the overall bandwidth, are already present in the 
smaller cluster. This is important since in the exact diagonalization study of (3xL)-ladders, 
we are restricted to relatively small (L<8) systems. 
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FIG. 3. Spinon excitation spectra for 7- and 8-site Heisenberg chains, (a) Single spinon excita- 
tion spectrum (S = 1/2) in a 7-site chain, (b) Two-spinon excitation spectrum (S = 1) in an 8-site 
chain, (c) Same as (a), but for a 19-site chain, (d) Same as (b), but for a 20-site chain. 
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For the (3 x 7)-ladder single spinon dispersion shown in Fig. 2(a), there appears an 
additional even-parity band, corresponding to the gapped spin degrees of freedom. It is 
separated from the odd-parity band by a gap ~ J, and its bandwidth is much smaller than 
that of the low-energy band. For the case of an even number of rungs (as shown in Fig. 2(b) 
for the case of (3 x 8)), the low-lying spin excitations now have even parity with respect 
to reflection about the center leg, and involve pairs of spinons with S=l and groundstate 
wavevectors k=± : | ± | = 0, ±7r. The finite size gap at k=±7r vanishes in the thermodynamic 
limit. Again, a comparison of the low-energy two-spinon band with that of the corresponding 
8-site single Heisenberg chain (Fig. 3(c)) shows qualitative (almost quantitative) agreement, 
supporting our conclusion that the low-energy degrees of freedom of the (3xL)-ladders can 
be mapped onto effective single chains. 

Finally, we analyze the finite size scaling behavior of the spinon velocity, v s = ^| k=0 , 
for the three- leg Heisenberg ladders (Fig. 2(c)). In the finite systems, we approximate 
the derivative by v s « ^Cgg/j^^TO (l even). Using L=2, 4, 6, and 8, we obtain with the 
approximate scaling form, v s (L) = v s (oo) + aL~ 2 + bL~ 4 , a bulk value of v s (oo) ~ 1.81 J, 
which is in excellent agreement (~ 1%) with that obtained by recent QMC calculations on 
clusters with up to 600 spins [TEfl . 

III. SINGLE HOLE IN A THREE-LEG LADDER 

We start the study of effects of doping by recalling the single electron non-interacting 
bandstructure. This takes the form of three overlapping bands (in the isotropic case t = t'). 
These can be classified according to their parity under the reflection operation (R) about 
the center leg. The two even-parity bands we denote as bonding (b) and anti-bonding (ab), 
and they have the dispersion relation 

e b ,ab( k ) = TV2t' -2tcosfc, (2) 
while the odd-parity band (or non-bonding band) has the form 
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£nb(k) = —2tcosk. 



(3) 



At half-filling, the chemical potential is /i — 0, and all three bands are partially filled. The 
Fermi surface consists of three pairs of Fermi points (±kp,A, A = ab, nb, b), arising from each 
band (or channel). 




anti-bonding (+) 
non-bonding (-) 
bonding (+) 



FIG. 4. Schematic plot of the non-interacting one-electron bandstructure of a (3 x L)-ladder. 
The parities of the bands with respect to reflection about the center leg are indicated by (±). The 
splitting between the bands is y/2t' . 

We consider the case of a single hole, and we start with a (3 x 7)-sample. The groundstate 
is now a total singlet (S=0) and non-degenerate. The wavevector dependence determines 
the single holon dispersion, again in a single chain effective model. In Fig. 5(a), the results 
are shown for a parameter value J/t = 0.5. The holon energy has absolute minima at k=±7r 
and a local minimum at k=0. This latter behavior is consistent with a finite size effect. For 
example in Fig. 5(c) we show the holon dispersion calculated in a 7-site single chain t — J 
model where a similar behavior is obtained. This in turn reflects the doubling of the period 
in the single Heisenberg chain dispersion, when the chain length goes to infinity. 

The three-leg ladder has reflection symmetry with respect to the center leg. The eigen- 
values of the corresponding operator, R, are ±1, i.e. even (odd) parity under reflection. The 
groundstate for the undoped ladder with (3xL) legs (L odd) has parity -1. The groundstate 
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FIG. 5. Hole excitation spectra for the isotropic (3 x 7)- and (3 x 8)-site t — J ladders, and for 

7- and 8-site t — J chains. In (a) and (b), the circles correspond to the even parity and the squares 
to the odd-parity channels respectively, (a) Single holon (8 = 1/21, S = 0) in a (3 x 7)-cluster. (b) 
Single hole excitation spectrum (S = 1/24, S = 1/2) in a (3 x 8)-cluster. (c) Filled circles: single 
holon excitation spectrum (5 = 1/21, S = 0) in a 7-site chain. Dashed line: same in a 19-site chain, 
(d) Filled circles: single hole excitation spectrum (S = 1/24, S = 1/2) in an 8-site chain. Dashed 
line: same in a 20-site chain. 



manifold of the same ladders with one hole has parity +1. Therefore we associate parities 
+ 1 and -1 with a single holon and spinon respectively. It follows that the parity of a single 
hole which is the product of those parities is -1, i.e. the hole goes into the band with odd 
parity with respect to R. This is the middle non-bonding band with odd parity. The in- 
terpretation is clear: the two bands with even parity in the t — J model combine to form a 
spin liquid which is insulating (ISL), and the initial doped holes go into the odd-parity band 
(or channel), and form a single Luttinger liquid (LL). So the non- interacting bandstructure 
Fermi surface is truncated from three sets of Fermi points to a single Luttinger liquid in the 
odd channel. Note that this Fermi surface truncation implies a form of spin pairing and a 
reduction of the low-energy spin degrees of freedom, but it does not imply the formation of 



10 



Cooper pairs. Note also that this Fermi surface truncation is not due to a breaking of trans- 
lational symmetry along the ladder since the spin-spin correlations in the ISL are purely 
short ranged. 
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FIG. 6. Spin-charge correlations for one hole in (3xL) isotropic t — J clusters with L = 4, 6, 8, 
and J/t = 0.5. The holon is fixed at the site indicated by the circle. 

White and Scalapino have performed DMRG calculations on large samples (3 x 16) with 



densities are well separated, indicating the decay of a single hole into a separated spinon 
and holon as expected for a Luttinger liquid. The Lanczos results for smaller clusters show 
similar results. In Fig. 6, the results for the spin-charge correlation are shown. If we add 
one hole to a (3xL) cluster with L even, then the groundstate has total spin S=l/2, and 
odd parity. Therefore if we take the state with S z = +1/2, we can calculate the local value 
of S z at each site for configurations with the holon fixed at the origin. In Fig. 6, we show 
these results for L=4, 6, and 8. The results show clearly that the spin is distributed over 
the whole cluster, and agree nicely with those of White and Scalapino |L(J . This behavior 
is fully consistent with a single Luttinger liquid interpretation. 

The Lanczos method allows one to examine also the low-lying excited states. This in turn 
raises the question of the energy gap between the Luttinger liquid with odd parity and the 



open boundary conditions [JTOf . They find upon doping a single hole that the spin and charge 
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TABLE I. Energies (in units of t) of the (3 x 6)-site isotropic t — J ladder with antiperiodic 
boundary conditions and J/t = 0.5. 
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energy 
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energy 





+1 
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vr/6 
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-10.24582 
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vr/3 
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vr/2 


+1 


-10.25131 
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+1 


-11.42219 
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+1 
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+1 


-10.25673 


2vr/3 


+1 


-11.49704 


7T 


+1 


-9.243230 
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+1 


-11.69586 
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-8.450508 


vr/6 


-1 


-9.987594 
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-11.59084 


vr/3 


-1 


-8.211651 


vr/2 


-1 


-10.40757 


vr/3 


-1 


-11.57428 


2vr/3 


-1 


-8.224553 
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-1 


-10.18534 


2vr/3 


-1 


-11.53402 


7T 


-1 


-8.403831 








7T 


-1 


-11.54535 



next channel with even parity. This latter channel should correspond to placing the hole in 
the spin liquid. From the experience with the two-leg ladder one expects that doping the spin 
liquid will lead to a Luther-Emery liquid with a bound spin-charge distribution for a single 
hole. Therefore if our expectation is correct, we should expect the lowest-lying excited states 
with even parity to show a very different spin-charge distribution. In Fig. 7, we show the 
corresponding distributions for the lowest eigenstates with even and odd parities of clusters 
with L=6. The eigenvalues are given in Table 1. The difference is clear. In the odd parity 
channel the energy dispersion is large, the minimum energy lies at k=±7r/2, corresponding 
to a half-filled Luttinger liquid. In the even parity channel, the energy dispersion is much 
less, and the minimum lies at values k=±57r/6 which we interpret as a hole entering the 
anti-bonding even parity band. As can be seen from Fig. 7, the instantaneous spin-charge 
distribution for this minimum energy even parity state is quite different, and shows a clear 
binding of the spin and charge. This behavior is fully consistent with that of a single hole 
in a Luther-Emery liquid and with the interpretation of a spin liquid in the even parity 
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channels. Also shown in Fig. 7 are the instantaneous spin-charge distributions for the set of 
lowest excited states with fixed wavevector and parity. These show intermediate behavior 
which we interpret as arising from (attractive) interactions between a spinon and a holon at 
wavevectors away from the Fermi surface. 
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FIG. 7. Spin-charge correlations for one hole in a (3x6)-site isotropic t — J cluster with 
J/t = 0.5. The holon is fixed at the site indicated by the circle. 



An important quantity is the single hole energy gap between the odd and even parity 
channels. The value we obtain in the (3x6)-cluster is »i 0.15t (or 0.30 J) - a value which is 
smaller than that of the spin gap in the ISL J). 

To summarize, the results for a single hole show that the minimum energy is in the 
non-bonding channel, and this channel forms a single Luttinger liquid. The even channels 
are gapped, forming an insulating spin liquid (ISL), and the minimum energy for a single 
hole in these channels lies higher in energy. 
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IV. TWO AND MORE HOLES 



We begin with the case of two holes. White and Scalapino have performed DMRG 
calculations for (3xL)-samples with a value of J/t = 0.5 [[ITJ. They find that the holes 
do not bind and in fact repel each other, so that the most likely configuration has the two 
holons as widely separated as possible, consistent with the open boundary conditions used 
in their study. This behavior can be immediately interpreted as that of two holes doped in a 
single Luttinger liquid (LL), while the remaining transverse channels form an ISL, consistent 
with the discussion given above. In view of the finite energy gap between the even and the 
odd parity states, there will be a finite density range within which this LL + ISL phase 
remains stable. We then conclude that the hole density in the White-Scalapino calculations 
(5 = 1/24) lies within this density range. 

The next issue is to determine the critical hole density, 5 C , that limits the stability of this 
LL + ISL phase. At first sight it would seem straightforward to use the single hole energy gap 
between odd and even parities to determine the critical value for the hole chemical potential, 
/i c , and thus 5 C : fi c = fi(S c ). However, we expect the even parity channels to evolve into a 
Luther- Emery liquid (LE) when doped, in analogy to the doped two-leg ladder. As a result, 
/i c will be determined by the energy of the two-hole bound state of the LE rather than the 
single-hole energy gap. 

In Fig. 8, we show the instantaneous hole-hole correlation function for various two-hole 
states of a (3x6)-cluster with APBC (J/t=0.5). These boundary conditions were chosen 
since they give a non-degenerate (i.e. closed shell) groundstate for non-interacting electrons 
||17|| . The groundstate has even parity and a total wavevector of k=0. The instantaneous 
hole-hole correlation function (Fig. 8(a)) shows the maximum weight at the largest rung- 
rung separation possible in this small cluster, but with both holes preferentially on the 
same edge. This behavior is similar to that found by White and Scalapino, and leads us to 
interpret this as a groundstate with LL + ISL character. Note that the hole density here 
(5 = 1/9), is larger than that considered by White and Scalapino. 
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FIG. 8. Instantaneous hole-hole correlations for two holes in a (3x6)-site isotropic t — J cluster 
with J/t = 0.5. The first hole is fixed at the origin (j=l), not shown in the figure, (a) groundstate, 
k=0. (b) excited state, k=7r. The sites on the first (outer) leg are labeled 1 < j < 6, on the second 
(center) leg 7 < j < 12, and on the third (outer) leg 13 < j < 18. 

There should however be an excited state of the cluster, also with even parity and S=0, 
which has the two holes in the LE. This state should have a total wavevector of k = tt, 
i.e. the same value as the undoped system jni. In fact, such a state appears in the cluster 
with an excitation energy of 0.1 lit above the groundstate. The instantaneous hole-hole 
correlations in this excited state (Fig. 8(b)) are also quite different from the groundstate, 
with a maximum for the next-nearest-neighbor separation, indicating that the two holes 
are bound in this state. Therefore we identify this state as that with two holes in a LE 
formed from the even parity transverse channels (i.e. bonding and anti-bonding bands). 
This identification is confirmed if we look at the average hole density on the central vs. the 
outer legs. The results quoted in Table 2 show a marked increase in the hole density on the 
central leg in the LE state. 

The binding energy of the two holes in the LE state can be estimated from the energy 
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TABLE II. Average hole densities on the central and the outer legs in the (3 x 6)-site isotropic 
t — J ladder with two holes and J/t = 0.5. 
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< n h > out 
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+1 


1.050053 
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tt/3 
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-1 
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vr/3 


-1 
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2vr/3 


-1 


0.978532 


1.042936 


7T 


-1 


0.884392 


1.231216 



difference to add the two holes together in the LE state (E LE ) and separately in the LE state 
(E LE ). This difference is an estimate of the binding energy, E b = 2E LE - E LE = -2.028t 
+ 2.453t = 0.4251 Note however that this value may well be a considerable overestimate 
because of finite size effects (see below). 

We now turn to the estimate of the chemical potential, In view of the small size 

of the clusters (i.e. mostly 3x6 sites) we need to take care about finite size corrections. In 
particular, in small clusters there is generally a large energy gain for total singlet states, and 
if one calculates \x through single hole additions one ends up subtracting the energies of states 
with total spin S=0 and S=l/2. In order to avoid this, we use two-hole additions since then 
the total spin quantum number remains unchanged; i.e. = |[E G (N h + 1) — E G (N h — 1)]. 
The result is shown in Fig. 9. The curve of /i(5) rises initially with 5, consistent with a 
repulsive interaction between holes in the LL + ISL phase. Our previous estimate of the 
excitation energy of the LE state with N h = 2 allows us to determine the chemical potential 
rise needed for holes to enter the LE phase at 0.055t. From Fig. 9 we estimate then for the 
critical hole density a value of 5 C ~ 0.13. Note that this value is calculated using a relatively 
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small cluster size, and it is not possible to estimate the finite size corrections. Therefore this 
value of 5 C is an estimate whose accuracy is hard to predict. 
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FIG. 9. Energy per site (circles) and chemical potential (squares) for a (3x6)-site isotropic t — J 
ladder with J/t=0.5. The critical hole density 5 C , beyond which the holes enter the LE phase is 
indicated by a dashed line. 

We conclude that the LL + ISL phase remains stable over a finite hole density range, 
0< 5 < 5 C , and in this range there is no hole pairing although a large fraction of the spin 
degrees of freedom are gapped in the ISL. Beyond this density range, 5 > 5 C , there is a LE 
channel in contact with the LL channel. The presence of holes in a LE channel leads to hole 
pairing and dominant superconducting or charge density wave correlations as in the case of 
the lightly doped two-leg ladder. 

In order to further explore the low-density LL + ISL phase, we calculated the single 
particle spectral function, A A (k, lo) = — ^Im[GA(k, u — ir/)], in the groundstate of the (3x6)- 
cluster with two holes. Here, A labels the linear combinations of c-operators, corresponding 
to the non-interacting shell structure shown in Fig. 5: Cj^a = \\ c j,i,a + V^Cj : 2,a + Cj,3,o], 
cj,nb,a = ^[cj,i l(T -Cj,3 iCT ], Cj^ a = \ [cj^„ - V2c j: 2,a + Cj,3,<x] ■ The results are displayed in Fig. 
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10 for (3x6) ladders with APBC and PBC. The energy region below the chemical potential 
corresponds to a transition to three-hole states. In a small cluster, such as the (3x6) under 
study here, the finite size effects are large. Thus to interpret the spectrum, one should keep 
in mind the non-interacting shell structure. We begin with the case of APBC. The two- 
hole LL + ISL state has then 16 electrons which in a non-interacting state would occupy the 
states fr(±f , ±f ), a6(±|), and n6(±|). Through the strong correlation interaction the state 
with electron pairs in the 6(±^) state will be admixed. In this way we can interpret the 
electron removal (or photoemission) part of the spectrum with the strong peaks to inserting 
holes into the 6(±f , ±^) and a&(±§ ) even parity channels. The state with an added hole at 
6(±|) is far removed from the Fermi level, and so it is strongly broadened. Turning to the 
nfr-channel, we see that the main weight is at (±|) as expected, but this is again broadened 
due to many-body effects. The resulting photoemission spectrum then shows signs of a 
Fermi surface in all three channels, and with kp-values in line with our expectations. 




0.0 1.0 2.0 0.0 1.0 2.0 

co/t Cfl/t 

FIG. 10. Spectral functions for a (3x6)-site isotropic t — J ladder with two holes and J/t=0.5. 
The poles have been given a finite width of i] = 0.02t, and the chemical potential is indicated by 
the dashed line. APBC (PBC) were used in the left (right) panel. 
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Clearly all three bands (b, ab, and nb) are partially occupied. This is illustrated by 
calculating the momentum distribution function, n^(k) = du>A\(k, to). The results for 
nA(k) are shown in Fig. 11. These show drops as k increases which one can interpret as 
giving estimates for k-p^x, but in a small system one cannot draw conclusions about the actual 
behavior for k ~ kp 5 A- The case of PBC has a different shell structure. The LL + ISL state is 
now a state with total momentum ±27r/3, with the 4 electrons in the n6-channel occupying 
k=0 (2 electrons) and k=±7r/3 (2 electrons) |l9j. The values of nA(k) for PBC in the LL + 
ISL state are also included in Fig. 11, and also confirm the interpretation of a partial filling 
of all three bands. 
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FIG. 11. Momentum distribution function for a (3x6)-site isotropic t — J ladder with two holes 
and «//i=0.5. Results with APBC and PBC were combined. 



However, at a true Fermi surface in a bulk system one must be able to add as well as 
to remove electrons at the chemical potential (or Fermi energy). The strong correlation 
condition clearly influences this part of the spectrum much more strongly, since as we pass 
to the undoped Mott insulator electron addition is totally forbidden (or allowed only on 
paying the Mott correlation energy gap), whereas in photoemission electrons can always be 
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removed. With this in mind we examine the electron addition part of the spectrum which 
corresponds to a transition to a one-hole state. From Fig. 10 we see that the main weight 
at low energy is in the n&(±|)-channel. Clearly this corresponds to a transition from the 
two-hole LL + ISL-groundstate to the one-hole LL + ISL-groundstate, and confirms that 
we have a Fermi surface in the odd-parity LL-channel. Next we look at the even-parity 
channels which form an ISL in the groundstate. Therefore we should expect that then 
electron addition in these channels at low energy will be forbidden. In fact, if we look at the 
6-channel, then this behavior is clear, and there is no Fermi surface in this channel. On the 
other hand, in the other even-parity channel, A = ab, we find two relatively strong peaks at 
k = ±|. These values actually correspond to a transition not to the minimum energy state 
of one hole in this channel but to an excited state. We interpret these peaks as sidebands, 
involving the creation of a S=l odd-parity magnon with k = or 7r in addition. Clearly it is 
too strong a statement to say that all processes which add an electron to a LL + ISL-state 
with even parity are forbidden, since sidebands in which an electron is added in the odd- 
parity LL channel and simultaneously odd-parity S=l magnons are created can have a total 
parity which is even and a total spin of S=l/2. We interpret the weight in the a6-channels 
then in terms of such processes which in turn may be enhanced by the finite size effects 
associated with these small clusters. The behavior of the PBC case is analogous. Here the 
main weight at low energy is the nb (k=7r/3) channel, corresponding also to a transition to 
the one-hole LL + ISL state. However, again there is weight in the antibonding channel 
which we interpret also as a higher energy sideband and a transition to a one-hole LE state. 
This is strongly admixed due to the enhanced stability of this state which has a filled shell 
character in the n6-channel as we discussed earlier. 

In conclusion, the single particle spectra in the LL + ISL state show a strong asymmetry 
between removing and adding electrons, and it is the latter process, which unfortunately 
is not easy to detect experimentally, where the effects of the Fermi surface truncation on 
approaching the Mott insulating state are most evident. 
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V. ANALYSIS BY MEAN FIELD THEORY 



In this section we analyze the properties of the three-leg ladder using the mean field 
description in the same spirit as previously done for the two-leg ladder [POR . For this purpose 
we introduce spinon and holon operators, / and b respectively, by replacing the electron 
creation and annihilation operators in the following way [PTj] . 

c i,u,a = fi,v,&bi,v and Q,!/,(t = b iu fi. U)t j (4) 

which lead to the local constraint J2a fl v vfi,v,o + &1iA,i/ = 1. The Hamiltonian can be 
reformulated in these operators, and can then be decoupled by introducing mean fields. The 
constraint is included by adding a term with Lagrange multipliers. Since this treatment is 
quite standard we do not go into details here. The system is not translationally invariant 
along the rung, so that the mean fields depend on the position of the bonds, resulting in 10 
independent mean fields: 6 hopping and 3 pairing mean fields, defined on nearest neighbor 
bonds (j, u;j',u'), 

Xj,v;j',u' = 2 J2a(fj,i/,crfj'y,u) 

Bj,v,j>y = (bj,vb], y ) (5) 

and a uniform Lagrange multiplier fi. We denote the mean fields along the outer two legs 
and the middle leg with the index 1 and 2, respectively, while the rung mean fields have 
the index 3. In order to avoid ambiguities we have to define the direction of the bond mean 
fields. We give here the convention using the example of x, 
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Xl — |Ect(/j>,(t/j+1,^) v — 1,3 

X2 = \ Y,a{f},2,afj + l,2,o) ( 6 ) 
X3 = ^Ea{fj,2,ak^) ^ = 1, 3 

The same convention is applied to the mean fields -B and A. The parity with respect to the 
reflection operator R allows us to separate the mean field Hamiltonian into the even (+) 
and odd (— ) parity channels, 

H mi = H? + + Hi f) + H<P 

+L[^ + i(2 x l + xl + 2A? + A|) + J'(xl + A|) (7) 

+2t{B lX i + B 2X2 ) + ^'Saxa] 
We define the corresponding combinations for the operators, operator combinations, 

/fc±,<r = yiz 52j(fj,i,o- ± fj,3,a)e lkr: > 



(8) 



with the momentum along the legs of the ladder. Then we can write the four terms of the 
Hamiltonian (in Nambu space for the spinons) as 
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H { _> = (-4*xi cos k - n)bl_b k - 



t /.t 



-4txi cos A; - fx -2v / 2t / X3 
— 2y / 2t'x3 — 4t%2 cos k — ji 



Jk+ 



\ b k,2 J 



H-^ — ^2a(fk-,cn f-k-,-a) 



-(2tB 1 + |Jxi) cos k — /i 



jJAi cos /c 



— I JAx cos /c 
(2tB 1 + |Jxi)cosA; + /i 



-^+ /) — (fk+pfk,2,pf-k+l,f-k,2,l) 



£*; A fe 
A fe -|_ fc 



/fc,2,T 

/-fc+4. 



where ^ and A fc are 2 x 2-matrices of the form 

-(2tB 1 + lJ X i)cosk- fi -2V2t'B 3 -^J' X 3 



-2^53 - ^J'xs -(2t5 2 + |J X2 ) cos k - n 



■fJAxcosA; -^J'A 3 



-3^J'A 3 _|ja 2 cosA; 



J k—,a 



f-k-,-0 j 



(9) 



(10) 



Note that we have neglected the terms with —riiUj/A in the exchange term, because in the 
mean field calculation they tend to artificially enhance the tendency towards pairing and 
favor a flux phase close to half-filling. 

We can now solve the single-particle problem of H mi and determine the mean fields self- 
consistently. For the groundstate we obtain the self-consistent equations by minimizing the 
groundstate energy of H m f in Eq.(4) with respect to all mean fields. In this description the 
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holons are Bose-condensed so that X^=i,2,3(frj,iA>) = 3$ where S is the doping concentration. 
In the following, we will use the same parameters as used in the numerical simulations, 
t = if = 2J = 21. 



A. Half filling 



For half-filling, 5 = 0, there are no holons and the Hamiltonian of the spinons is invariant 
under the following SU(2) transformation 



(f! 

-r 



( 



u 



\ 



JjM,a 
\ fj,v,-tr ] 



( 



\ 



f'J 

ft 



(11) 



which reflects the constraint at half-filling. The absence of an up-spin corresponds to the 



presence of a down-spin and vice versa [22]. Note that due to the SU(2) invariance of the 
Hamiltonian the self-consistent solution of the mean field is not unique, but U corresponds 
to a rotation in the mean field space {xi,2,3, A.1,2,3} which leaves the groundstate and the 
excitation spectrum unchanged. 

In this formulation various properties of the three-leg ladder can be interpreted at least 
on a qualitative level. The Hamiltonian yields six distinct bands in Nambu space for the 
spinons (Fig. 12 ). For the symmetric channel the spinons have a gapped spectrum where 
the lower two bands are completely filled. On the other hand, the antisymmetric channel has 
gapless excitations which give a spectrum analogous to that of a single chain. It is now easy 
to compare the spin excitation spectrum of the mean field and the exact diagonalization. 
For the mean field case the excitation corresponds to a particle-hole excitation of the spinon 
gas. The gapless spectrum is that of a single Heisenberg chain and is exactly what we obtain 
for the same type of mean field treatment of a single chain. The other bands are gapped 
and describe a spin liquid state coexisting with the gapless system. The agreement with the 
numerical simulation is qualitatively very good. 
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FIG. 12. Spinon bands at half filling obtained by the mean field theory. All states below the 
chemical potential (located at ek = 0) are occupied by a spin-up and a spin-down spinon. The 
gapped bands correspond to the even-parity states, while the gapless bands in the center are in 
the odd-parity channel. 

B. A single hole 

We now insert a single hole into the half-filled system, i.e. we remove one spinon and add 
one holon. This immediately breaks the SU(2) symmetry, because the groundstate energy 
of the holon has to be minimized which yields a constraint on the hopping mean fields Xi- 
Three holon bands appear, two even parity bands and one odd parity band, with energies 
given by 

e k+ = -2t(xi + X2) cos k - fi 

± v /4t 2 (xi-X 2 ) 2 cos 2 A; + 8t 2 xi ( 12 ) 
e fc _ = —4:txi cos k — ji. 
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Obviously the lowest holon state can be found in the lower of the two even-parity bands. 
Thus, the single hole groundstate can be obtained by removing a spinon at the Fermi level 
of the odd-parity band and by inserting a holon at bottom of the lower even-parity holon 
band. This results in a change of the parity of the groundstate compared with that of the 
half-filled case, in agreement with the numerical result. 




FIG. 13. Holon bands for a single hole obtained by the mean field theory. The outer 
bands are the even-parity bands (bonding and anti-bonding) and the center band has odd par- 
ity (non-bonding). 



The single particle spectrum is incoherent, because the excitations are composed of the 
annihilation of a spinon and the creation of a holon. We define the single hole Green's 
function as a matrix 



G ul/ >(k,u) 



dte-™\T(cl(t)c ky , a (0)) 



(13) 



where the indices v and v' denote the three legs of the ladder (T the time-ordering operator). 
Decomposing the c-operators into spinon and holon parts we obtain the convolution 



GW(*,<") = \ll\ da/Gls(k + q, U + <^G b w ,(q,u>') 



(14) 
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with 




FIG. 14. Low-energy spectrum for a single hole obtained by the mean field theory. 
The spectrum is obtained by diagonalizing the Green's function with respect to the 
leg indices. We show in Fig. 14 the lowest one-hole states for given momenta and parity 
analogous to the numerical results in Fig. 5. The spectrum changing the parity with respect 
to the groundstate of the half-filled system has gapless excitations at momentum k — 
and 7r. This corresponds to the LL component of the three leg ladder. The states which 
keep the parity of the groundstate, however, are gapped with minima at k = and n. 
In our mean field calculation the latter excitations lie above the LL spectrum, in contrast 
to the numerical result. The reason is that in the mean field treatment these excitations 
are almost exclusively carried by the spinon part, i.e. we remove a spinon in the lower 
even-parity spinon bands, while the holon remains in the lowest even-parity band. In the 
numerical calculation this spectrum lies clearly lower and cannot be identified directly with 
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spinon excitations which are higher in energy. Therefore, we conclude that our mean field 
calculation overestimates the splitting of the holon bands. Nevertheless, we can interpret 
both types of excitations consistently with the numerical calculations. In particular, we 
emphasize that the gapless spinons are in both treatments in odd-parity channel giving a 
consistent picture of the symmetries of the ground state and the low-lying excitations. 



The analysis of our numerical data led us to the conclusion that there would be a finite 
range of doping close to half-filling where the LL state would coexist with the ISL. Beyond a 
critical doping S c the ISL would be doped and the whole system would turn into a LE where 
all channels open a spin gap. We might hope to describe this behavior within our mean field 
picture. However, in this treatment the mean field LL state is immediately unstable against 
the formation of a LE state, and there is no obvious transition between two different phases. 

As soon as we introduce a finite concentration of holes in the system, the holons are 
described by a Bose-condensate in the mean field groundstate. Consequently, the single- 
particle Green's function consists of a coherent and an incoherent part due to the presence 
of the condensate, 



where C vv > is a constant proportional to the doping concentration 8. Thus, we find conven- 
tional Fermi liquid quasiparticles in the single particle spectrum. This has implications on 
the pairing in the doped ladder. The exchange term in the t-J model allows to scatter spin- 
singlet pairs of electrons between the even-parity and odd-parity channel. The relevant terms 
are derived from the exchange terms along the two outer legs, JXa(Sj,i • + &i,3 ' Sj+i^), 
and become in electron operator formulation, 



C. Finite doping 



G^iKu) = C uu ,Gl u ,(k,Lu) + G2 c ,(k,u;) 



(16) 



H. 



+- — 



J2k,k',q £<r=U v k,k',q [ c k+q+,a C k'+,-a C k-,-(TCk'+q-,a- 



(17) 



a Ck+,-aCk'+q+,a\- 
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with Vk,k', q = — J(cosg + |cos(A; — k')). The even-parity channel corresponds to the ISL 
phase which is not populated by holes at small doping. However, via virtual scattering of 
hole pairs into the even-parity channel an effective pairing interaction for the holes in the 
odd-parity channel is generated, in lowest order given by 



V k , k ' = -2i Efc" / duv Kk i^ k+k »G s ++ {k" ', uj)G f ++ (-k", - 



^)Vk',k",k'+k" 



-UJ] 



= — [ duj(9 cos k cos k' cos 2 k" + sin k sin k' sin 2 k")G+ + (k" , u)G+ + (— k",—c 
2 k „ J 

(18) 

where G+ + (k, lu) is the spinon Green's function in the even-parity channel corresponding to 
G+ + = | (Gii + G{ 3 + G31 + G33) from Eq.(15). In weak coupling theory this leads to the 
instability equation (ladder approximation) in the odd-parity channel, 

1 = -k B Tj2Vk,kG—(k,u n )G^(-k,-u n ) (19) 

k,n 

for finite temperature (G = G\\ — G13 — G31 + G33, and u n are the fermionic Matsubara 

frequencies). Since the Green's function has a standard quasiparticle pole as shown above in 
Eq.(16), we find a usual Cooper instability and a non- vanishing transition temperature for 
any finite coupling and density 5 (note that G is proportional to 5). Thus, holes in the odd- 
parity channel would be paired in the groundstate. This is in contrast to the expectation 



for a LL where a critical coupling strength exists below which the LL remains stable p3 
Thus, this instability of the LL for any non- zero concentration of holes is a deficiency of the 
mean field description. On the other hand, it clearly indicates the trend towards pairing 
due to the coupling of the LL to the insulating or doped spin liquid which occupies a part of 
the spectrum. The RVB correlation of the ISL provides a pool of "preformed pairs" which 
triggers the instability of the LL towards the LE |24f . 

Although there is no obvious (non- zero) critical concentration where the character of the 
groundstate changes within the mean field treatment, we observe an interesting crossover 
between two regimes in the BCS mean field for the holes, 
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5 

FIG. 15. The superconducting mean fields for finite doping. The solid lines correspond to 
the mean field calculation. The dashed lines are fits of the crossover assuming an underlying 
y/5 — <5 c -dependence of A' with S c « 0.034. In the lower figure a linear fit is possible for A' 2 . 

= (Ci^) = (blb])(frf n ) = BijAij. (20) 

The behavior of A' and A' 2 as a function of 5 is shown in Fig. 15. The crossover indicates 
a critical concentration 5 C pa 0.034 (for t = 2 J) found by extrapolation of the square root 
dependence of A' (oc y/5 — S c ). For < 5 <C 5 C A' is rather small as can be understood 
from Eq.(15) and (16). For small doping the dominant contribution to Vk,k' originates from 
k" m ±7r/2 the momenta for the spin gap in the even-parity channel and the Fermi points 
in for the odd-parity holes are close to ±n/2 as well. Therefore the cosine-term in Vk,k> 
contributes only little. This changes with increased doping where both the Fermi points 
and the lowest even-parity spinon energies are located at momenta gradually shifting away 
from ±7r/2. This leads to an effective enhancement of the pairing interaction as an effect 
of doping which yields the pronounced crossover whose symmetry aspects we discuss below. 
We expect that the same tendency is at work in the real system where in the very small 
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doping regime the effective interaction is sufficiently reduced to avoid pairing in the LL. Only 
with increased doping both the enhanced attractive interaction and the density of states at 
the Fermi level drive the system into the paired state. 

Finally we would like to characterize this crossover by considering the symmetry of the 
gap function in momentum space. Thus, we introduce transverse momenta k± along the 
rung. The bonding, antibonding and non-bonding states as given in Eq.(2) and (3) can be 
related to momenta k± by using the following form for the dispersion, 

efc,fc ± = —2t(cosk + cos/c_i_), (21) 

resulting from fixed boundary conditions along the rung. With Eq.(2) and (3) we find the 
correspondence: b — > k± = ±7r/4, nb — > k± = ±tc/2 and ab — > ±3ir/4. Then we obtain for 
the momentum dependent gap function, 

(22) 

= -(2A' 1 sin 2 (A;_ L ) + A' 2 sin 2 (2A;j.)) cos k - A' 3 sin(2A;_ L )(sin(A;_ L ) + sin(3A;j.)). 

We can use this form now to show the phase of the gap function in the first Brillouin zone 
(BZ). We observe that the gap function is basically positive along the /c-direction (\k±\ <C it) 
and negative along the transverse momentum direction k± (\k\ <C n). This is essentially the 
structure of a "d^^-wave" pairing function. This is also reflected in real space where the 
gap function along the legs is positive (A'^A^) and along the rungs is negative (A' 3 ) as 
shown in Fig. 15. It is interesting to compare the position of the nodes in the two regimes 
separated by S c . For < 5 <C 5 C we find that the gap functions have practically the same 
magnitude on all bonds. With this property Eq.(22) leads to nodes which correspond exactly 
to the [110]-direction in the BZ and we might consider this state as "purely" d x 2_ y 2-wave like 
(see Fig. 16). On the other hand, for the regime 8 > 5 C the gap functions have the relation 
A^ < A' 2 = — A' 3 which yields nodes clearly shifted away from the [110]-direction. We may 
consider this as an admixture of an (extended) s-wave-like component to the d-wave gap, 
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although, strictly speaking, the underlying symmetries are not present in the ladder system 
to justify the distinction between s- and d-wave pairing. 
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FIG. 16. Nodelines for the BCS gap, Eq. (22), in the first Brillouin zone. The dashed line 
connecting the empty circles shows the nodelines for 5 = 0.01. The dotted-dashed line connecting 
the x-marks gives the nodelines for S = 0.06. We define the transverse momenta k± by b — > ±7r/4, 
nb — >■ ±7r/2, and ab — >■ ±3ir/4. 

The difference of the two regimes may therefore be interpreted in the following way. The 
mean field superconducting state in the small doping regime is mainly carried by the even- 
parity channel and the odd-parity channel participates weak through proximity. Contrary 
to the exact numerical diagonalization, in the mean field treatment we cannot avoid the 
population of the ISL by holes strictly even at very small doping concentrations. The three- 
leg structure is apparently wide enough to create a LE state with a gap structure which 
is approximatively d-wave like. In the larger doping regime, however, the LL of the odd- 
parity channel acquires its own hole pairing so that an additional symmetry lowering occurs 
reflecting the ID nature of the LL and generates an s-wave like additional component as seen 
in the shift of the nodelines (Fig. 16). Since there is no real symmetry distinguishing s- and d- 
wave from each other this transition appears only as a crossover when the interaction among 
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holes become strong enough in the (odd-parity) LL. This behavior is again in qualitative 
agreement with the interpretation obtained from our numerical calculation. 

VI. CONCLUSIONS 

The three-leg ladder is specially interesting because in a certain sense it combines two 
quite different elements - an odd-parity channel which behaves like a single chain, and two 
even-parity channels which behave similar to the two-leg ladder. This result is already 
foreshadowed in the undoped limit where the low-energy degrees of freedom can be mapped 
onto a single Heisenberg S=l/2 chain with longer range interactions, and the remaining 
transverse spin degrees of freedom have a substantial energy gap. In this work, we show 
that the initially doped holes enter only this odd-parity channel and form a Luttinger liquid, 
so that the system has two distinct components - the conducting Luttinger liquid in the 
odd-parity channel, and the insulating spin liquid in the even-parity channel. This phase 
we denoted as LL + ISL, and the numerical evidence from the Lanczos diagonalization of 
small clusters reported here and the DMRG calculations of White and Scalapino [[KJ , that 
such a phase exists up to a critical hole density, 8 C , is we believe quite clear. In addition, we 
introduced a mean field approximation scheme which gave similar but not identical results. 
Initially, holes enter only the odd-parity channel which is gapless in the undoped system. 
However in the mean field approximation, a small gap develops in the odd-parity channel 
upon doping. This - as we discussed - reflects the inadequacy of the mean field description 
of the Luttinger liquid. 

The LL + ISL phase has unusual properties. First of all, we note that the different parity 
channels of the original Fermi surface behave quite differently, so that only in the odd-parity 
channel is there a Fermi surface. The truncation of the Fermi surface in the partially occupied 
even-parity channels is not a consequence of a breaking of translational symmetry since the 
spin order here is purely short range. Rather it is a consequence of the proximity to the Mott 
insulating phase which in this channel is ISL. This truncation of some partially occupied 
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bands is a clear violation of Luttinger's theorem. Usually if one approaches a Mott insulator 
which has AF order, then one may proceed through incommensurately ordered phases which 
progressively truncate the Fermi surface, but do not violate Luttinger's theorem. However 
this option is not available if one approaches an RVB Mott insulator which is an ISL. However 
the example of the three-leg ladder shows us that here also a partial truncation of the Fermi 
surface is possible, but now it violates the Luttinger theorem. 

This LL + ISL phase has certain features in common with a recent proposal by Geshken- 
bein, Ioffe, and Larkin (GIL) for the underdoped spin gap normal phase of the cuprates. 
They argue that the spin pairing implied by the spin gap did not cause immediately hole 
pairing. Instead they broke up the Fermi surface into two distinct parts - a fermionic part 
and a paired bosonic part. The latter they argued should have infinite mass to prevent 
conductivity from these bosons, and they associated this with the van-Hove singularity. In 
our specific example, the spin pairs are also insulating, but the origin lies in the proximity to 
the Mott insulating phase. In both models, the coexistence of fermionic and bosonic degrees 
of freedom leads to processes whereby a fermionic Cooper pair can scatter in and out of the 
bosonic channels, which however lie at higher energy. In a LL, there is no Cooper instability 
for an infinitesimal attraction, and a finite attraction is required for pairing. For this reason 
we believe a LL + ISL phase is possible in the three-leg ladder. By contrast, as we discussed 
above, the mean field description of the odd-parity channel has a true Cooper instability 
due to the holons being Bose condensed, and as a result hole pairing occurs at arbitrarily 
small hole densities. In the case of two dimensions, where the Fermi surface channels or 
patches near to the saddle points, (±7r, 0) and (0, ±7r), first become paired and insulating, 
it is clearly crucial whether the remaining fermionic part of the Fermi surface has a Cooper 
instability (as assumed by GIL) or not. A Cooper pairing instability leads to hole pairing 
in the groundstate as in our mean field description. 

The experimental examination of the Fermi surface evolution has been made by ARPES 
( Angle- Resolved-Photoemission-Spectroscopy). However, as we have seen, the inverse pro- 
cess would be more illuminating, but it is much more difficult to realize experimentally. In 
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the electron addition process, as we have shown, the Fermi surface truncation is very evi- 
dent, although sidebands do allow an electron addition on the even-parity channel which is 
in the ISL state. Nonetheless the total weight at low energies will vanish as the hole doping 
vanishes, and will be concentrated mainly in the fermionic part of the Fermi surface. 

In the three- leg ladder at dopings beyond <5 C , the ISL is converted into a doped Luther- 
Emery liquid, and a LL + LE phase occurs. This phase will show hole pairing and power 
law correlations in the CDW and singlet superconductivity channels. As the mean field 
approximation shows, the pairing is in an essentially d-wave channel. In the mean field 
theory, the quantum phase transition at 5 C appears as a crossover where the hole pairing 
increases rapidly as the hole density increases. 
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